Computing Congruence Primes Between a 
Newform With Integer Coefficients and the Old 

Space 



Abstract 

In this article, we put together some known theoretical results and 
the fact that certain computations can be done efficiently in SAGE 
to come up with a fast algorithm for calculating congruence primes 
linking a newform with integer coefficients (i.e. a newform associated 
to an elliptic curve) with the old space at the same level. 



1 Introduction 

Let be a positive integer. Throughout this article, let / be a newform of 
weight 2 and level N with integer coefficients. Let X be some subspace of 
S'2(ro(A^)). Recall that we call a prime p a congruence prime between / and 
X if there exists some cusp form g in X with integer Fourier coefficients such 
that ttnif) = o-nig) mod p for all n, (in which case, we write f = g mod p). 
Congruence primes play an important role in number theory. For instance, 
they were used extensively in the work of Ribet |Rib90j and Wiles [Wil95j . 
among others, leading to the proof of Fermat's Last Theorem. 

Let Mx be the set of forms in X with integer coefficients. Suppose 
that one can efficiently compute a finite set {gi . . . gr} of forms in Mx that 
generate Mx as a Z module. (This is possible, for example, when X is the old 
space.) In this article, we give a method for quickly computing congruence 
primes between / and such a space X. 

We became interested in computing congruence primes when trying to 
test two conjectures of Agashe regarding cancellations in the conjectural 
Birch and Swinnerton-Dyer formula for the elliptic curve associated to / and 
their relation to congruences between / and old forms, which we describe 
in detail in Section 5. Agashe was able to test these conjectures using the 



SAGE S"'"ll| command congruence_number() for curves of conductor up to 



1000, at which point computations became too slow to gather more data. 
The method we present is fast up to level 1500 and is fast up to much higher 
levels when the levels are nonsmooth (e.g. if A/^ is a product of two primes). 
In Sections 2 and 3, we present the algorithm for the general subspace X. 
In Section 4, we specialize to the case when X is the old space, where it is 
especially easy to compute an integral spanning set for Mx- In Section 5, 
we explore the algorithm's motivation and main application and compare it 
to congruence_number(). 

2 Deciding if a Given Prime is a Congruence Prime 

Recall that / is a newform of level N with integer coefficients, X is a subspace 
of S'2(ro(A^)), and Mx is the set of forms in X with integer coefficients. In 
this section, we give a criterion for deciding if a given prime p is a congruence 
prime between a newform / and X. 
Let 

[SL2{Z) : ro(iV)] [SL^jZ) : ro(jV)] - 1 
[ 6 N J- I • J 

The following is a consequence of Theorem 1 in |Stu87| . 

Theorem 1 (Sturm). Let h and g he two cusp forms in S'2(ro(A^)) and let 
p he a prime. If an{h) = an{g) mod p for all n up to B then an{h) = an{g) 
mod p for all n. 

The integer B above is called the Sturm hound of ro(A^). If is a cusp 
form, then let v{g) denote the row vector whose components are the first B 
coefficients of g. Let {gi, . . . ,gr} be a set of cusp forms with integer Fourier 
coefficients whose Z-span is Mx- Let M be the integer matrix whose i-th 
row is v{gi) for i=l, . . . , r. 

Lemma 1. The newform f with integral coefficients is congruent modulo p 
to a form in X if and only if the v{ f) is in the row space of M modulo p. 

Proof. If v{f) is in the row space of M modulo p, then we must have a row 
w = Civ{gi) in the row space of M such that Wn = v{f)n mod p for all 

n up to B. Thus an(X]i=i ^i9i) = (^n{f) mod p for all n up to B. Theorem 
[1] applies so that an(Yl\=i ^i9i) = '^n{f) mod p for all n. 

Conversely, if there exists a form g = Y^i=i ^i9i ™ such that an{g) = 
ttnif), for all n, then v{g)n = v{f)n for all n up to B. Then v{f) = 
Yll=i '^i'^idi) mod p, so that v{f) is already in the mod-p span of M . □ 
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In section m we describe how one can efficiently compute M associated 
to X in tlie case that X is the old space, so that Lemma [1] can be used to 
test if a prime is a congruence prime linking a newform to the old space. 

3 Finding all Congruence Primes 

In this section, we describe an algorithm for finding all congruence primes 
between the newform / and the subspace X of 5'2(ro(A^)) when one is given 
a finite set of forms that generate Mx as a Z module. 

Let E be the elliptic curve over Q associated to /. Consider the map 
(j) : Xq{N) — > Jo{N) given by (/> : P i— )• (P) — (oo). The map (p induces a 
surjective morphism 

Xo{N) ^ Jo{N) ^ E. 

The degree of this composite morphism is called the modular degree of E. 

The following is theorem 3.11 in |ARS06) . 

Theorem 2 (Agashe, Ribet, Stein). Suppose f is a newform of level N with 
associated elliptic curve E, and f is congruent modulo p to a cusp form in 
its orthogonal complement in 52(ro(A^)). Then either p divides the modular 
degree of E or p'^\N . 

The modular degree can be computed very efficiently. For a history of 
computational approaches of calculating the modular degree, see the intro- 
duction in |Wat02| . In SAGE, Watkins's algorithm [Wat02] for computing 
the modular degree of an elliptic curve is called by modular_degree(). 

Lemma [T] provides an efficient test to check if a prime is a congruence 
prime between / and the oldspace, while Theorem [2] gives a finite list of 
candidate primes that contains all congruence primes. Together, we have 

Algorithm 1. Given a newform f of level N and a spanning set for X of 
forms with integer coefficients, this algorithm computes all congruence primes 
between f and X. 

1. Compute the matrix M as described in Section\^ 

2. Compute the modular degree and generate finite list of primes {pi-.-Pn} 
such that Pi divides the modular degree or pf\N . (Recall that by The- 
orem [3 these are the only primes that could be congruence primes 
between f and X .) 

3. For each pi, check if the v{f) is in the row space of M modulo pi. 
4- Output the list of primes for which the answer in step three is yes. 
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4 Computing the Matrix M when X Is the Old 
Space 



Let X be the old space at level A^. We need to find a set of forms {gi, . . . , g^} 
with integral coefficients that generate Mx over Z. Recall that for some d\N ^ 
the d-degeneracy map is given by 

oo oo 

Pd:Y.<f)'i'^Y.^i^f)l''- (4-1) 

1=1 i=l 

For each prime p\N , all old forms originating from any level dividing ^ 
factor through S'2(ro(^)) and pass to 5*2 (Fq (A'")) by the 1 and p degeneracy 
maps. Thus, to compute an integral spanning set for the old space at level 
N , it suffices to compute an integral basis for the spaces 5*2 (Fo(^)) for p\N 
and pass these forms up to level by the 1 and p-degeneracy maps. This 
discussion leads to the following algorithm. 

Algorithm 2. Given a level N , this algorithm returns an integer matrix M 
whose rows are v{gi) where {gi, . . . ^g^} is a set of cuspforms with integer 
coefficients that spans the old suhspace of S2{Tq{N)). 

1. Let M be an empty matrix. 

2. For each prime p\N : 

(a) Compute an integral basis of S2{To{^)) : 

In SAGE, the command integral _basis() can be used to compute 
an integral basis for 5*2 (Fo(-^)). 

(b) Compute the 1-degeneracy images from spaces S'2(Fo(^)).' 

Let g be a form in a basis obtained in step (a). By equation \4-l\ 
Pi{g)=g. So for each form g in each the basis obtained in step 
(a), augment to M the vector v{g), which consists of the first B 
coefficients of g. In SAGE, the command coefficients () can be 
used to compute the Fourier coefficients of a cusp form. 

(c) Compute the p-degeneracy images from spaces S'2(Fo('^)).' 

By equation \4.1\ the vector v{/3p{g)) is a vector of length B whose 

i-th component is zero if p\ i and v{g) i_ if p\i. For each g in the 

p 

basis obtained in step (a), augment to M the vectors v{(3p{g)) . 

3. Output the matrix M . 
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The reason that Algorithm 2 is fast in practice is because the command 
integral_basis() in step 2(a) is very efficient. In Section 5, we indicate why 
Algorithms 1 and 2 compute congruences between / and X faster in SAGE 
than congruence_number() when X is the old space. 



5 Motivation and Application 

In this section, we describe the problems that us to be interested in com- 
puting congruences between newforms associated to elliptic curves and the 
old space, and we indicate why our method is fast compared to congru- 
ence_number() for this case. Let E/Qhe a, modular elliptic curve associated 
to the newform / of level N. The right hand side of equation conjectural 
Birch and Swinnerton Dyer formula, 

L^{E/Q,i) _ \mEm\ ■ R{E/q) ■ ripCp 



r! n \E, 



tors 



has been studied extensively (see the introduction in |Lor| ). In particular, 
in regard to the relationship between the Tamagawa product IIpCp and the 
order of -EtorsCQ), M. Emerton |Eme03) showed that when is prime, IIpCp = 

|^tors(Q)|. 

When A^ is not prime, |-E'tors(Q)| need not equal IIpCp. In fact, this first 
occurs when A^ = 42. D. Lorenzini [Lor] has shown that if £ > 3 is a prime 
number, then the order of the i!-primary part of -E'tors(Q) divides IIpCp. 

In the other direction, Agashe |Aga| has conjectured that: 

Conjecture 1 (Agashe). If an odd prime i divides UpCp, then either i divides 
the order of -Etors (Q) or the newform f is congruent modulo i to a form in 
the old space. 

The following partial result towards Conjecture [1] is given in |Aga| . 

Proposition 1. Letl be an odd prime such that either i f A^ or for all primes 
r that divide A^, ^ f (r — 1). If i divides the order of the geometric component 
group of E at p for some prime p||A^, then either E[i] is reducible or the 
newform f is congruent to a newform of level dividig N/p (for all Fourier 
coefficients whose indices are coprime to N£ ) modulo a prime ideal over i in 
a number field containing the Fourier coefficients of both newforms. 



In fact, in |Aga the author conjectures more specifically that: 
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Conjecture 2 (Agashe). // an odd prime i divides Cp for some prime p, 
then either i divides the order of E'torslQ) or the newform f is congruent 
modulo i to a form in the subspace generated by degeneracy map images of 
newforms of levels dividing N/p. 

We became interested in computing congruence primes when trying to 
test Conjectures 1 and 2. Agashe was able to test these conjectures using 
the SAGE |S+ll| command congruence_number() for curves of conductor 
up to 1000, at which point computations became too slow to gather more 
data. The method we present is fast up to level 1500 and is fast up to much 
higher levels when the levels N are nonsmooth (e.g. if A'^ is a product of two 
primes) . 

The algorithm used in congruence_number() in SAGE is designed to 
calculate the congruence number between any two disjoint subspaces X and 
Y of 5'2(ro(A^)). It does so by computing the integral structures, Xjs and 
Yjs, of the modular symbol spaces associated to X and Y, and finding the 
index of Xjs + Yjs in its saturation - a process that requires: 

1. computing a Z-basis for the modular symbol spaces corresponding to 
X and Y 

2. calculating the Hermite normal form H of the matrix associated to 
Xis + Yis 

3. finding the Hermite normal form M of 

4. taking the determinant of M. 

This method can be modified to check only if a prime p is a congruence 
prime by performing steps 2-4 modulo p, but we do not know whether or not 
such an approach would detect all congruence primes without false negatives, 
and we do not know how efficient the modified algorithm would be. 

In the special case that X is the span of a newform with integer coeffi- 
cients and Y is the old space, the method we described for computing congru- 
ence primes is faster than congruence_number() for two reasons. First, the 
results in Section [2] allow us to perform operations modulo certain primes. 
Second, as described in Section [U when Y is the old space, the problem of 
finding an integral basis for Y can be reduced to computing integral bases 
of the cuspidal spaces at lower levels, which is already efficient in SAGE. 
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